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Abstract 

The family of harmonic Hilbert spaces is a natural enlargement 
of those classical L^-Sobolev space on which consist of continuous 
functions. In the present paper we demonstrate that the use of basic 
results from the theory of Wiener amalgam spaces allows to establish 
fundamental properties of harmonic Hilbert spaces even if they are 
defined over an arbitrary locally compact abelian group G- Even for 
Q = this new approach improves previously known results. In this 
paper we present results on minimal norm interpolators over lattices 
and show that the infinite minimal norm interpolations are the limits 
of finite minimal norm interpolations. In addition, the new approach 
paves the way for the study of stability problems and error analysis for 
norm interpolations in harmonic Hilbert and Banach spaces on locally 
compact abelian groups. 
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1 Introduction 



Babuska ^ introduced the concept of periodic Hilbert spaces in order to 
study universally optimal quadrature formulas. Subsequently Prager (16] 
has studied in detail the relationship between optimal approximation of linear 
functionals on periodic Hilbert spaces and minimal norm interpolation. In a 
more recent paper Delvos |6j has introduced the notion of harmonic Hilbert 
spaces for the real line M and discussed the interpolation problems over an 
infinite uniform lattice for £^-data in M. 

In the present paper we introduce the concept of harmonic Hilbert spaces 
over a locally compact abelian (LCA) group Q, establish the properties of 
these spaces, and provide a characterization of minimal norm interpolators. 
The use of Wiener Amalgam spaces provides a clear understanding of the 
subject. This paper opens new viewpoints for the study of optimal norm in- 
terpolation in various types of harmonic Banach spaces over locally compact 
abelian groups. 

The article is organized as follows. In Section |21 we introduce the notation 
and state basic facts of harmonic analysis on compact abelian groups that we 
essentially need for proving the main statements. The next section gives a 
short overview on Wiener Amalgam spaces and contains a summary of some 
of their important properties. In Section H] we introduce weight functions 
and derive the reproducing kernel of harmonic Hilbert spaces. Then we show 
in Section El how the concept of Riesz-basis leads to certain principal shift 
invariant subspaces. In particular, it follows that the bi-orthogonal system 
of the Riesz sequence of uniform translates of the kernel is generated by the 
same family of translates of the so-called Lagrange interpolator. Section [HI 
contains the core results of the paper, namely the explicit derivation of mini- 
mal norm interpolations. In the final section we discuss the approximation of 
the minimal norm interpolation by optimal interpolations of finite samples. 

2 Preliminaries and Notation 

Let ^ be a locally compact abelian group and Q its dual group with the nor- 
malized Haar measures dx and d'j, respectively. The dual group Q is defined 
as the set of all characters j : Q ^ T, i.e., continuous homomorphisms from 
Q into the torus T with pointwise multiplication as group operation, and the 
standard (compact-open) topology. 
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For 1 < p < oo, L^-spaces are denoted by their usual symbols. For 
/ G L^{Q), the Fourier transform JF is defined as 

^/(7) = /(7)= / f{x)^dx. 
Jg 

There are many excellent textbooks on abstract harmonic analysis such as 
[T7j and 113] that comprise a rigorous treatise of Fourier transform on LCA 
groups. Here, we only state those facts that are most useful for the develop- 
ment of the presented results. 

Using the Pontryagin duality theorem which allows to write (7, x) or (x, 7) 
unambiguously for 7(x), the Fourier inversion formula for h = f E L^{Q) 
writes as 

f{x)=J'-'h{x)= / /(7)(7,x)ci7- 
Jg 

The range of the Fourier transform, which is an injective and bounded linear 
mapping from to C^, the space of continuous functions that vanish at 
infinity, is denoted by JFL^. The convolution theorem states that 

Hf*9)=Hf)H9), f.9eL\ 

where " *" denotes the usual convolution. By transport of the norm ||JF/||jfli = 
11/11^1, the normed space [TL^ ^ || • turns into a Banach algebra with 

respect to pointwise multiplication. Due to Plancherel's theorem one has 
TV- = * L^, which ensures that TL^ is a dense subspace of C° endowed 
with the sup- norm || ■ ||oo. 

Whenever a Banach space (X, || ■ is continuously embedded into an- 
other Banach space (y, || ■ i.e., there exists a constant C > such that 
ll^^lly ^ for all x G X, we write X '-^Y . 

Throughout this paper we use the symbol for the translation operator 

Txfiy) = f{y - x), x,yeg. 

For convenience, all index sets, sums, and lattices on Q and Q used in the 
sequel are assumed to be countable. 

Let A be a subgroup of Q. According to [Hj, we call A a lattice if 
the quotient group Q/A is compact. The lattice size s(A) is defined as the 
measure of a fundamental domain of A in Q, i.e., we choose a measurable 
set ?7 C ^ such that every x E Q can be uniquely written as x = A -|- m for 
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some A G A and u & U, cf. ^3]. Then the lattice size s(A) is the measure 
of U. For an equivalence class in ^/A we write x = {x + A}. Let dx denote 
the normalized Haar measure on Q/A. For / G L^{Q), we then have Weil's 
formula 

(^J2f{x + X))dx. (1) 

AeA 

The annihilator of A is the subgroup 

A^ = {xeg\xW = i, AgA} 

and Weil's formula becomes in this context the identity 

/ = s{A^) [ ( J2 /(7 + X))di . (2) 

Note that the dual group A of A is naturally isomorphic to Q/A-^, [T7] . 
For (ca) G ^^(A), the sum XIaga '^^(7, A) is a Fourier series on Q/A-^ and 
Plancherel's theorem ^\ yields 

/ |$^c,(7,A)'d7 = $^|cAr- (3) 

•^e/A^ ' AeA A6A 

For Q = R and A = aZ, a > 0, this is just the classical Fourier series 
expansion of periodic functions. In this case we simply have A-*- = 

We make use of these fundamental results in the proof of Theorem 15.21 



f{x)dx = s(A) 



g/A 



3 Wiener Amalgam Spaces 

Wiener Amalgam spaces on LCA groups and their properties under Fourier 
transform are studied in a series of papers starting with jH]. We refer to 
[TH 113] for a compact survey on their properties. Here, we only provide a 
short overview over the facts needed in the present context. 

First we introduce the continuous version of Wiener Amalgam spaces. 
Let ip he a. non-zero, non-negative function on Q with compact support, i.e., 
supp('?/') C Q where Q C G is compact. Given any Banach space {B, \\ ■ \\b) of 
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functions on Q on which JFL^ acts by pointwise multiphcation, i.e., !FL^-B C 
-B, the Wiener Amalgam space W{B, L'^) is defined by 

W{B, L^) = {fe Bloc I ||/lk(B,L.) := ( ^ 11/ • T.^lllrfx)'^' < oo } (4) 

for 1 < g < oo (with obvious modifications for q = oo). The expression 
/ G Bloc means that the function / can at least locally be measured by the 
_B-norm hence the integrand in the definition of || • ||vi/(_b,L9) is a well-defined 
non-negative function. As stated in [TT] , these spaces are Banach spaces and 
do not depend on ip (up to equivalence of norms). 

An equivalent, but discrete, definition of Wiener Amalgam spaces uses a 
so called jFL^-bounded, uniform partition of unity (for short BUPU), that is 
a sequence of non-negative functions {ipi} corresponding to a sequence {gt} 
in Q such that 

1. supp(V^i) C gi + Q, 

2. sup, ^{j\{gi + n)n {gj + r^) } < oo, 

3. EMa)^^ 

4. ipi are bounded in JFL^, i.e., H-i/^jUjFLi < C < oo for all i. 

With the help of such a BUPU we define the discrete Wiener Amalgam space 

WiBJ'^) = {fe Bloc I := ( Yl IIMIIb)'^' < oo } (5) 

i 

for 1 < g < OO. Note that again, W{B,i'^) is a Banach space independent 
from the partition of unity (up to equivalence of norms). In all cases of this 
paper we will use a lattice A in ^ instead of an arbitrary sequence {gi}- 

In contrast to the L^-spaces, where, for instance, on M'^, inclusion fails 
for different p, Wiener Amalgam spaces enjoy the so-called coordinate-wise 
inclusion, i.e., if Bi^oc ^ B2,ioc and qi < ga, then W{BiJ''^) ^ W{B2,£''^). 

Wiener Amalgam spaces can also be defined with additional weight func- 
tions that are described in the following section, cf. 
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4 Reproducing Kernel Hilbert Spaces 



We first state sufficient conditions on weight functions that are used to define 
harmonic Hilbert spaces. It is beyond the scope of this work to elaborate on 
these conditions. 

A strictly positive and continuous function w is called a submultiplicative 
(or Beurling) weight function on Q if 

^(71 + 72) < w(7l)^^^(72) for all 71,726^, (6) 

It satisfies the Beurling- Domar non-quasianalyticity condition 17, VI, sect. 3] 
if 

00 

(BD) J2^~'^^'^Sw{n'y) < 00 for all 76^. (7) 

n=l 

The standard example of such weight functions are weights of polynomial 
type on M'^ such as 

u;,(7) = (l + |7|)^^(l + |7p)^/2^ 7GR^ 

for s > 0, or subexponential weights such as 

w(7) = e"l^l'(l + |7|)% 7GM^ 

for a > 0, < 5 < 1, and s G M. For detailed studies on such weight 
functions on the Euclidian space M'^ we refer to [71 CHI • 

Throughout this paper, we assume that G L'^{Q) which implies w~'^ G 
W{C^,£^). Obviously, this is satisfied for the above example Wg whenever 
s > d/2. In contrast to [Uj where the author included the box function as a 
possible weight leading to band-limited functions that we discuss separatly in 
Example 2, we restrict our discussion to submultiplicative weight functions 
satisfying (BD) in order to obtain stronger results and general statements 
valid for arbitrary lattices. From our point of view the band-limited case 
should be seen as a limiting case, requiring sometimes separate arguments. 

Among others, the submultiplicativity of w in conjunction with the as- 
sumption G L^{Q), allows to apply the convolution theorem for Wiener 
Amalgam spaces jHJE], in order to derive the following crucial property: for 
any lattice A there exist positive constants a and b such that 

a<Yl ^~'h + x)<b, leg. (8) 
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The upper bound follows from the fact that the Haar measure for A-*-, i.e. 
/i = X^xga^ belongs to W{M, where M denotes the space of bounded 
measures, while on the other hand the submultiplicativity implies that w~'^ G 
W{C'^,f ). Since the periodization of w"^ equals /i * w"^ C W{M,i°°) * 
C W{C,i^) = C\G) the upper bound is valid, see [11]. Since 
is strictly positive the periodicity implies that there is a strictly positive 
lower bound. 

We define the harmonic Hilhert space corresponding to the weight w as 

'H^{Q)=r-'Ll{G) (9) 

with inner product 

Jg 

Lemma 4.1. 'Hw{G) is a Hilhert space. 

Proof. Left to the reader. □ 

In the case of the weight functions w^, Tiw^ is known as the Sobolev space 
of fractional order s 

Remark. It has been shown in ^ that if G is subadditive, i.e., 
w(7i + 72) < C{w{-ii) + w{-i2)) , 7i, 72 e ^ , 

then is a Banach convolution algebra. As a consequence, Tiw turns into 
a Banach algebra with respect to pointwise multiplication. This is the case 
for Ws, s > d/2, and subexponential weights. 

We now state a result which is fundamental for sampling functions in Hw 
along some lattice A. 

Theorem 4.2. If w'^ G L'^{G), then 

Proof By virtue of Corollary 7 of [HI, we have J^-^Ll = W^T-^Ll^f). 
We emphasize that this result is a consequence of the submultiplicativity 
and (BD). Now, since G L^, an easy application of the Cauchy-Schwartz 
inequality implies ^ L} which in turn results in the following inclusions 
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(by the Riemann-Lebesgue Lemma). Hence, by the coordinate- wise inclusion 
properties of Wiener Amalgam spaces we obtain 



□ 

As a consequence, TCw{G) is a reproducing kernel Hilbert space (RKHS). 
This is a Hilbert space consisting of continuous functions in which the point 
evaluation functionals are continuous. Therefore, by virtue of the Riesz Rep- 
resentation theorem, for each x & Q there exists a unique function, say fc^;, 
in the RKHS, such that the point evaluation at x of any function / in the 
RKHS can be performed by means of the inner product with k^. In the case 
of Tiyj, that is 

f{x) = {f,k^)^, f eHyj{Q) ,x eQ . 

The kernel k{x, y) = kx{y) defines a continuous function on ^ x ^ containing 
all the informations about the scalar product. A first detailed survey on 
RKHS goes back to 1.2.^. 

In the present situation the kernel of Tiw turns out to consist of transla- 
tions of a single function. 

Proposition 4.3. TiwiQ) is a reproducing kernel Hilbert space with kernel 
k{x, y) = (j){x — y), where 

That is, f{x) = (/, Tx(f))w for all f G Hw and x E Q . 

Proof. Define by = w~^. Since G L'^, (p G Tiw ■ Next we compute 

{f.TA)n, = [jh)f^i7Wi7)d7 
Jg 

= / hl){^^l)w~'^ilWil)d'y (10) 
Jg 

f{l){x,l)d-i = /(x), 

the last step following from the Fourier inversion theorem (note that / G L^). 
The uniqueness of the kernel completes our proof of the fact k^ = Tx4>. □ 
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Remark. The fact that Tiw C° can also be derived immediately from 
(imp . The use of Wiener Amalgam spaces, however, reveals the important 
property that the sequence of samples of any function in Tiui with respect to 
any lattice A is square summable. This property has led to the definition of 
so-called puzzles in ITS]. 

Corollary 4.4. For any lattice A the mapping 

Q:f^ (/(A))AeA (11) 

is bounded from TiwiG) to £^(A). 

Proof. Choose a fixed lattice A' of Q with fundamental domain Q! . Then 
the family {Ta/XhIa'cA' where x denotes the characteristic function, forms a 
BUPU. We now consider the norm of W{C,l'^) with respect to this BUPU. 
For any A' G A', there is at most a finite number of points of A in \' + Vt' , 
say nxi- These numbers are uniformly bounded. Therefore, ||(/(A))||£2 < 
sup(?7.A')||/||vy(c/2). Hence, the result follows from Theorem 14.21 □ 

We will later see that the mapping Q is surjective. 

5 Riesz Basis 

Riesz bases are a well-established concept in Hilbert space theory 0. 

Definition 5.1. A family of vectors {hn} in a Hilbert space Ti is called a 
Riesz sequence if there exist bounds < a < b < oo such that 

a||c||^2 < II ^Cn/in||^ < H^Wp (12) 
n 

for all sequences c = (c„) G i"^. 

Riesz sequences generalize the concept of orthogonal sequences as one 
can see from the following properties. Let us call V to be the closed linear 
span of {hn}. For every Riesz sequence {hn} there exists a unique dual or 
bi-orthogonal sequence {hn} in V such that 

{hn, hm) = Snm, Ti, m G N , 
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and the orthogonal projection Py from Ti. onto V is given by 

Pyh = ^^(/i, hn)hn = ^^(^5 ^n)^n , for all h E H , (13) 



cf. A Riesz sequence {/?.„} obviously constitutes a (Riesz) basis for V. 

The following statement is a standard result in Fourier analysis, cf. [TH] . 
For the sake of completeness we include the prove. 

Theorem 5.2. For any lattice A 0/ ^, the sequence {Tx0}agA ^-5 Riesz 
sequence in L'^{Q) if and only if there exist positive constants a, h such that 



a<Yl I<^(7 + X)l'<& 



a. e. 



(14) 



Proof. We compute 



I 2 Planchorcl 
I2 ~ 



AeA 



AeA 

L |<^(7)n X^Ca (A, 7)1^^7 

f|<^(7 + x)ri$^CA(A,7 + x)| 



s(A 
s(A 
s(A 



xeA^ 



AeA 



S/A-L 



xeAJ 



AeA 



g/AJ 



X] |0(7 + X)n X^ca(A,7)|^ ^7 



xeAJ 



AeA 



□ 



•^^/^^ xeA^ AeA 

Since this holds for all (c)a G ^^(A), the statement follows from (0). 

A similar result holds for Tiw^Q) instead of L?'{Q). 

Theorem 5.3. For any lattice A of Q, the sequence {TA^jAeA "is a Riesz 
sequence in TiwiQ) if o.'^^d only if there exist positive constants a, b such that 



xeA-L 



a. e. 



(15) 
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Proof. Analogue to the proof of Theorem 15.21 □ 

RecaUing Condition (jH)) on the weight, we finally obtain that {Ta^Iasa 
forms a Riesz basis for its closed linear span, the so-called spline-type space 

\4(0) = span{TA0|AG A} = { / = caTa0 | (ca) G £'(A) } 

AeA 

which is a closed subspace of TiwiQ)- 

We now look for the dual basis of {Ta^Iasa in Tiw- Let us define 

= (16) 

Note that the denominator is a A-'--periodic square integrable function. Since 
any such function has a unique Fourier series presentation with ^^-coefficients, 
we can easily see that if) belongs to Va(0). Similar to the proof of Theorem 
15.21 we compute 

(Ta0,Tv^)^ = 5a,a', A,A'gA. (17) 

Hence, {TaV^Iasa is the dual basis of {Ta^Iasa in Va(0). It follows from (jTH]) 
that every function in Va(0) can be written as 

/ = Ta^) Ta^ = ^ /(A) TaV^ . (18) 

AeA AeA 

As a consequence, every function in Va(0) is completely determined by its 
samples on A. 

Because of the reproducing kernel property of 0, (fT7|) is equivalent to 
V^(A - A') = 5a,a' , A,A'G A. 

Therefore, if) is also called the Lagrange interpolator for A. 

It is important to note that the bi-orthogonal system of {TA</)}AeA is again 
generated by the same translates of a single function, namely the Lagrange 
interpolator. 

Remark. Following Theorem 15.21 Va(0) is also a closed subspace of L^. 
Similar to above, the dual basis of {Tx^'jAeA in L"^ is generated by the dual 
atom ip2 given by 



^2(7) 



E,eAx|0(7 + x)l^ 
11 



Example 1. In ^21; the authors give a detailed study of interpolation and 
stable reconstruction of functions in the Sobolev space Hw^iM.'^) from samples 
taken over a lattice of varying lattice size. 

Example 2. A well-known example is the space of band-limited functions 

B={fe L'iR) I supp(/) C [-1/2, 1/2] } 

endowed with the L^-inner product. It might be seen as a Harmonic Hilbert 
space for the weight function 



WiXj 




■1/2,1/2], 
-1/2,1/2] 



although the weight function does not satisfy the conditions (jHI) and ((Tj). 
Nevertheless, results similar to those above can be obtained. For instance, 
5 is a RKHS with the so-called sinc-kernel 

sin Tcx 

smc(x) = 

nx 

which is the inverse Fourier transform of the box function 

1 XG [-1/2,1/2], 



[-1/2, 1/2]. 



Note that ^ = The set of integer shifts of the sinc-kernel constitutes an 
orthonormal sequence. It is even an orthonormal basis of B. Consider the 
lattice A = aZ with a < 1. The corresponding annihilator is A-*- = a~^Z, 
and it can easily be seen that 



fcgZ 



is not bounded away from zero. In particular {Taksmc}k£Z is an overcomplete 
basis system, a so-called frame jB] and not a Riesz basis, in accordance to 
Theorem 15.21 The fact that the periodization that is not bounded away 
from zero almost everywhere, leads to a frame, is a general property and 
characterization of such frames as shown in j3] . 
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6 Minimal Norm Interpolation 

The problem of minimal norm interpolation in a harmonic Hilbert space on 
the real line has been discussed by Delvos, 0. In the present section we 
study the corresponding problem in a more general setting by integrating 
results stated above. 

Theorem 6.1. Let A be a lattice in Q and (ca) G £^(A). The interpolation 
problem /(A) = C\ for all X & A has a unique minimal norm solution in 
'Hw{Q). It coincides with the interpolating element in Va(0) with = 

Proof. Since Theorem 14.21 shows that T-Cw{Q) "-^ W{C,£'^){Q), the mapping 

/ G HUG) - (/(A)) G f{A) 

is well-defined and bounded, cf. Corollary 14.41 Due to Theorem 15. 3( {Tx0} 
is a Riesz basis for Va(0) which, by duality, consists of all functions of the 
form 

AeA 

for some (ca) G and the Lagrange interpolator ip which is the dual Riesz 
atom for {Ta0}. Due to the Lagrange property of ip it follows that 

/(A)=ca, AG a, 

for any such / G V\{(f)). Hence the mapping 

/ G - (/(A)) G f{A) 

is bijective. In particular, there exists a unique element fc G Va(</>) with 
/(A) = Ca. Since Va(0) is closed, we can split Tiw^Q) into the direct sum 

n^{G) = v^{<p)®vt{<P), (19) 

where = { / e | /(A) = 0, A G A } . Assume that some g G 

Hw{G) interpolates c on A. Then g — fc E V^(0) and we obtain 

Ml = \\{9- fc) + fcWi = \\9- fcWi + ml > wfcWi ■ 

Hence, fc is the unique minimal norm element in TiwiG)- D 
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Corollary 6.2. The minimal norm interpolation of the sequence (/(A)) for 
some f G 'Hw{Q) sampled on a lattice K in Q coincides with the orthogonal 
projection of f onto Va(0), say Pf. 

Proof. This follows immediately from 

P/(A) = (T,</., Pf) = {PTx<P, f) = {T,<j>, f) = /(A) . 

□ 

7 Orthogonal Projection 

In the final section we show that the minimal norm interpolation for a finite 
number of lattice elements converges to the minimal norm interpolation for 
A when increasing the number of lattice elements. 

Let A be a lattice of Q. We denote by {Ap} a nested sequence of finite 
subsets of A with [JpAp = A. We define 

Vp{4>) = span{ Tx(f)\XeAp}. 

We recall that | Ta0 | A G Ai? } is obviously a Riesz basis for Vp{(f)) whose 
dual basis is given by { PpTxip | A G A^r } where Pp denotes the orthogonal 
projection onto Vf(0). It is obvious that PpTxip = for all X ^ F. 
For a fixed element c = (ca) G £^(A) we have seen that 

9 = Y1 ^^^^^ 
AeA 

is the minimal norm interpolation on A. Set Cp = {cx)x^\p. By the same 
arguments used in Theorem 16.11 we easily deduce that the minimal norm 
element for in HwiG) is given by 

X£Ap XeAp AeA 

In other words, the minimal norm interpolation for cp is just the orthogonal 
projection of the minimal norm interpolation g for c. 

Exploiting the nested structure of the subspaces, a standard argument in 
wavelet theory, e.g., [201; immediately implies 

\\g — gpllw ^0 for F increasing . (20) 

Since TiwiQ) is continuously embedded in C(^), the convergence in holds 
true also for the sup-norm. 
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